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We show that there are two reasons why the partial width for the transition 
ri(T(3S') — > 7 + xw(l-P)) is suppressed. First, the spin-averaged matrix ele- 
ment I(3S\r\lPj) is small, being equal to 0.023 GeV -1 in our relativistic calcu- 
lations. Secondly, the spin-orbit splittings produce relatively large contributions, 
*^H' giving I(3S\r\lP2) = 0.066 GeV -1 , while due to a large cancellation the matrix 



element I(3S\r\lPi) = —0.020 GeV 1 is small and negative; at the same time the 



magnitude of I(3S\r\lPo) = —0.063 GeV 1 is relatively large. These matrix ele- 
> ■ 

ments give rise to the following partial widths: T2(T(35) — > 7 + Xfc2(l-P)) = 212 eV, 



ro(T(3S') — > 7 + Xbo(l-f)) = 54 eV, which are in good agreement with the CLEO 



m 

On 
O 

and BaBar data, and also to ri(T(3S') — > 7 + Xfci(l-P)) = 13 eV, which satisfies the 
O ■ 
<N ; BaBar limit, ri(exp.) < 22 eV. 
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.?_' I. INTRODUCTION 

In recent years, several new bottomonium states were discovered due to studies of radia- 
tive decays [1-4]. In [1] CLEO has observed the T(1D) in the four-photon decay cascade, 
T(3S)^ 1 + Xb (2P), X6 (2P)-> 7 + T(l£>), T(LD) -> l + Xb(lP), Xb(lP) -> 7 + T(lS% 
and later this state was observed by BaBar in another four-photon cascade via the T(2S) 
[2]. In 2008 a new state, rjb(lP), was discovered by BaBar, first in radiative decay 
T(3S) — > 7 + r]b(lS) [3] and then in T(25) — >■ 7 + rjb(lS) [4]; later 776(15*) was confirmed 
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by CLEO [5]. Moreover, new or more precise data on different radiative transitions, like 
T(3S) -). 7 + Xb(n 3 Pj) (n = 1,2), X6 (1P, 2P) -> 7 + T(15), and Xfc (2P) -+ 7 + T(25), 
were presented in Refs. [6-9]. 

This new experimental information is of a special importance for the theory to provide a 
better understanding of the role of relativistic and spin-dependent effects in bottomonium, 
and may be used as a test of different models and approximations. There are a large number 
of papers devoted to radiative decays in bottomonium [10-14], and a comparison of different 
results was already presented in [12-14], where the predicted partial widths are shown to be 
rather close to each other for most radiative El transitions and to agree with the existing 
experimental data. The only exception is the radiative decays T(3S) — ¥ 7 + Xb(lPj) {J — 
0, 1, 2), which are discussed in detail in [14]. Their partial widths are defined by the matrix 
element (m.e.) I(3S\r\lPj) = (T(3S)\r\l 3 Pj) (J = 0, 1,2) and below we shall also use the 



spin-averaged m.e., denoted as I(3S\r\lP). 

These m.e. strongly differ in the nonrelativistic (NR) and relativistic cases, even within 
the same model. The predicted transition rate ri(T(3S f ) — > 7 + Xfei(l-P)) varies in a wide 
range, (3 — 110) eV [14] and is in many cases larger than the experimental width: T\ = 
(33 ± 10) eV from the CLEO data [8]); a smaller value r\ = (10lg) eV was measured by 
BaBar [9]. Moreover, even in the models which predict a small partial width Ti, their other 
two rates, Tj(J = 0,2), do not agree with the experimental values [15]. Therefore, the 
ratio of the transition rates, nn = tw-y-L en 1+x ^ 5^ 1 must be considered an important 
characteristic, which is small in experiments: r^o ~ 0.5 from the CLEO [8] and r^o ~ 0.2 
from the BaBar data [9] . 

The m.e. I{3S\r\lPj) may differ several times in NR and relativistic calculations, even 
within the same model or while different static potentials are used [11, 16]. In Rcf. [17] 
the suppression of this m.e. was shown to be quite strong in the NR limit for the power- 
law potentials V(r) ~ r a with — 1 < a < 2. Since in bottomonium, even for T(3S), the 

2 

relativistic corrections are not large, -^ < 0.1, one may assume that this fact occurs because 
of the different asymptotics of the wave functions (w.f.) of the Schrodinger and relativistic 
equations. 

An interesting result was obtained in Ref. [16], where for the NR Hamiltonian the partial 
width T 2 = r(T(3S f ) — > 7 + Xb2^P)) decreases ten times, if instead of the Cornell potential 
with a (static) = constant, the Wisconsin potential which takes into account the asymptotic 



freedom behavior of the vector strong coupling, is used. This result reminds of the situation 
with the dielectron widths of T(nS) (n = 1,2,3), where agreement with experiment is 
reached only for the potential with the asymptotic freedom behavior of the strong coupling 
[181. 



However, even for this kind of potentials the spin-averaged m.e. I(3S\r\lP) appears to 
depend on the freezing (critical) value of the vector strong coupling used. In this paper we 
consider gluon-exchange (GE) potentials with two different values of a CT - lt . 



It is also evident that since the m.e. I(3S\r\lPj) is small, it may strongly depend on 
other small effects, in particular, on the spin-orbit interaction used. Here we show that due 
to the spin-orbit splittings the m.e. I(3S\r\lPj) acquire corrections of the same order as 



the value of the spin-averaged m.e. I(3S\r\lP), and a large cancellation takes place in the 
m.e. with J = 1. Here in our calculations we use the relativistic string Hamiltonian (RSH) 
[19], which was already tested in a number of papers, devoted to different bottomonium 
properties [20]. 

II. RADIATIVE DECAYS 

Electric dipole transitions between an initial state % = 3 3 S\, and a final state / = 1 3 Pj, 
are defined by the partial width [10-14], 



r(i^ 1 + f) = -ae 2 Q E*(2J'+l)S i r f \£ if \ 2 , (1) 

where J' — Jf, I' — If, and the statistical factor S^j- = S^ is given by 



Sj} = max(J,J')< J l J ' \ (2) 

1 s l J 

and for the transitions between the n 3 Si and m 3 Pj states with the same spin 5 = 1 this 
coefficient Sfj = 1/9. 

The RSH is simplified in the case of bottomonium, where in the Hamiltonian the string 
and self-energy corrections can be neglected because they are very small, < 1 MeV. Then 
the original form of the RSH with the static potential 

V B (r) = ar- (3) 

3 r 



H= P_±JJh +C0 + y Bir y (4) 
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Here nib is the 6-quark pole mass, while the value of u is determined from the extremum 
condition |^ = 0, which gives u = \Jp 2 + m 2 , being equal to the kinetic energy of a 6 
quark. Substituting this u into Eq. (4) one arrives at the spinless Salpeter equation (SSE): 



H = 2y/p 2 + m 2 b + V B (r), (5) 

The kinetic term occurring in (5) is widely used in relativistic potential models [21-23], 
however, as compared to constituent potential models, the RSH has several important dif- 
ferences. 

1. By derivation, the mass of the b quark in the kinetic term cannot be chosen arbitrarily: 
it must be equal to the pole mass of a b quark, which takes into account perturbative 
in a s (mb) corrections. In two-loop approximations m^pole) = m&(m&)(l + 0.09 + 
0.05) [24], where the second and third numbers come from the a s and a 2 corrections, 
respectively. In our calculations m^pole) = 4.83 GeV is used, which corresponds to 
the conventional current mass fhb{fhb) — 4.24 GeV. 

2. H , as well as the mass M(nl), does not contain an overall additive (fitting) constant. 

3. The string tension a = 0.18 GeV 2 , used in the RSH, cannot be considered a fitting 
parameter, because it is fixed by the slope of the Regge trajectories for light mesons. 

4. In the GE potential the asymptotic freedom behavior of the vector strong coupling 
o:b(V) is taken into account, being expressed via the "vector" QCD constant Ab, which 
is not a fitting parameter but defined by the conventional Aj^g according to the rela- 
tion: A B (ra/ = 3) = 1.4753 A ] ^(n / = 3) and A B (n f = 5) = 1.3656 A ] ^(n / = 5) [25]. 
On the other hand, the value of A-jjg{nf = 5) is fixed by the known value of a s (Mz) at 
the scale Mz = 91.19 GeV. Here a s (Mz) = 0.1191 is used, which in two-loop approx- 
imation gives Ajjg(nf = 5) = 240 MeV and correspondingly, Ab(ti/ = 5) ~ 330 MeV. 

Thus our scheme of calculations appears to be very restrictive in the case of bottomonium 
and only small variations of the fundamental parameters are admissible. However, some 
uncertainty comes from the value of the freezing constant, a B {r — > oo) = a cr n, which 



properties are discussed in Ref. [26]. Here we use the vector coupling in the range 0.49 < 
c^crit < 0.60. Then for a given multiplet nl the centroid mass M cog (nl) coincides with the 
eigenvalue M{nl) of the SSE: 



2 \ /p 2 + m 2 b + V B (r) 



ip ni = M(nl)if n i. (6) 



For this relativistic equation the NR limit and the so-called einbein approximation may 
also be used and in both approximations a good description of the bottomonium spectrum 
is obtained, even for the higher states [20]. For most radiative decays (in bottomonium) 
the m.e. like I(mS\r\nPj) and I(nPj\r\mS) differ only by 10 — 20% in the NR and rel- 
ativistic cases, with the exception of the transitions T(3S') — > 7 + Xbj(l-P)- bn this case 



our calculations give I(3S\r\lP) = 0.007 GeV l in the NR case, being ~ 3 times smaller 



than I(3S\r\lP) = 0.023 GeV" 1 for the SSE (here a crit = 0.49 was used). Notice that for 
a stronger GE potential with a crit = 0.60 these spin-averaged m.e. appear to be larger: 



I(3S\r\lP) = 0.011 GeV" 1 in the NR case and 0.036 GeV" 1 for the SSE. 

Since the same static potential is used for the SSE as in the NR case, such a difference 
between the m.e. may be explained by two factors: the different asymptotic behavior of the 
w.f. of the SSE and Schrodinger equations, and also a smaller value of the w.f. at the origin 
for the Schrodinger equation as compared to that for the SSE. However, it is known that the 
w.f. R n s(r), as well as the derivative R' 1P (r) for the IP state, diverge near the origin for the 
SSE (these divergences are discussed in details in Ref. [16]) and the calculated values of the 
w.f. (or its derivative) at the origin are obtained with the use of a regularization procedure. 
This regularization introduces a theoretical error, which is estimated to be < 10%. 

In Table I we give the m.e. I(3S\r\lPj), calculated here for the SSE and in the NR limit, 
together with their values from second paper of Ref. [13] (table 4.16) for the NR and the 
relativistic variant RA, where a scalar confining potential, as in our calculations, is used. 

Comparison of the m.e. presented in Table I shows that 

1. In Ref. [13] for the relativistic variant RA the m.e. I(3S\r\lPj) is < 4 times larger 
than in the NR case; a similar result is obtained here for the spin-averaged m.e., where 



I(3S\r\lP) = 0.023 GeV" 1 for the SSE and is an ~ 3 times smaller value 0.007 GeV" 1 
in its NR limit. 



2. Corrections SI so (J) = I(3S\r\lPj) — I(3S\r\lP), due to the spin-orbit potential, have 



TABLE I: The m.e. I(3S\r\lPj) (in GeV ) in the relativistic and NR cases 



Transition 


NR 


RA a ) 


NR 6 ) 


SSE 




[13] 


[13] 


this paper 


this paper 


<3S|r|LP 2 ) 


0.016 


0.063 


0.047 


0.066 


(3S|r|lPi) 


0.011 


0.063 


-0.033 


-0.020 


(35|r|lP ) 


0.004 


0.063 


-0.073 


-0.063 



a > Given numbers refer to the variant RA [13], where a scalar linear potential is used. 
' Given numbers refer to the NR limit of the SSE Eq. (6) with the same potential Vb(v) and 

a CTit = 0.49. 



TABLE II: 


The partial widths 


r(T(35) -> 


7 + XbO-Pj)) (m eV) 




Transition 


i 


RA 


NR 


SSE 


exp. 


exp. 




(MeV) 


[13] 


this paper 


this paper 


CLEO [6] 


BaBar [8] 


r 2 (T(3S)^7 + X6(lP 2 )) 


433.5 


195 


108 


213 


157 ±30 


216 ±25 


r 1 (T(3S)^ 7 + Xb (lPi)) 


452.1 


134 


36 


13 


33 ±10 


< 22 


r (r(3S)^7 + x6(LPo)) 


483.9 


54 


72 


54 


61 ±23 


55 ±10 



a relatively large value, e.g. 5I so (J = 2) = 0.043 GeV , being almost two times larger 



than I(3S\r\lP) in the spin-averaged case (see Eq. (9) below). 

3. In Ref. [13] the splittings between the m.e. I(3S\r\lPj) with different J are much 
smaller than in our calculations. 

4. In the spin-orbit potential we take the strong coupling a so (fi) = 0.38, which is close to 
the value a so (//(2P)) used for the Xbj(2P) states (this value was extracted in Ref. [23] 
from the experimental masses of the members of the Xbj{2P) multiplet). Our calcula- 
tions here show that the nondiagonal m.e., like (nP\V so \mP) (n ^ m,n = 1,2,3), are 
of the same order or have even larger values than the diagonal m.e. (2P|V^ (r)|2P). 

The calculated El transition rates are presented in Table II together with their values 
from Ref. [13]; they correspond to the m.e. from Table I. 

In the relativistic case our transition rates appear to be very close to those from the 
BaBar data [9]. Even in the NR case, due to large spin-orbit corrections, the calculated 
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partial widths do not contradict the CLEO data [8]. 

We make some remarks on the contribution 5I so to the m.e. I(3S\r\lPj) from the 
spin-orbit potential, V so (r) — L ■ S V so (r), for which the splittings a so (nP|lP) = 
(nP\V so \lP), (n = 2,3) are taken as for the one-gluon exchange interaction, i.e., neglecting 
the second order corrections in a s (fi) (it may be shown that the second order corrections 
are negative and small, ~ —0.7 MeV). In this approximation we find 

a so (nP, IP) = 7-^{4a so (r _3 ) nP> ip - a(r' l ) nPtlP }, (7) 

where we take a so = 0.38, which provided a good description of the fine-structure splittings 
for the Xb{2Pj) multiplet. To determine the corrections to the w.f. of the Xbj(lP) states, 
the potential V^ is considered as a perturbation and the following mass differences between 
the centroid masses are used: 

M cog (2P) - M cog (lP) = 360 MeV, M cog (3P) - M cog (lP) = 640 MeV. (8) 

Notice that the correction from the 3P state is not small, while the value of the centroid 
mass M(3P), M(xb(3P)) ~ 10.54 GeV, is taken from the recent ATLAS experiment [27]. 

For the SSE, the splittings a so (2P, IP) = 12 MeV and a so (3P, IP) = 10.2 MeV, were 
calculated and in the NR limit their values are ~ 10% smaller. Then the nondiagonal m.e. 
I(3S\r\lPj) with the "spin-orbit" corrections can be presented (in GeV -1 ) as 



I(3S\r\lPj) = I(3S\r\lP) + SI so (J) : 



8l so (J) = 0.033 O I(3S\r\2P) + 0.016 O I(3S\r\3P), (9) 



where O = -2, -1, +1 for J = 0, 1, 2 and I(3S\r\lP) = 0.023 GeV -1 for the SSE (relativistic 
case) and 0.007 GeV -1 in the NR limit. To obtain the m.e. presented in Table I, we use 



also the spin-averaged nondiagonal m.e.: I(3S\r\2P) = —2.54 GeV 1 and I(3S\r\3P) = 
2.64 GeV -1 . 

III. CONCLUSIONS 

For the El radiative transitions, T(3S*) — > 7 + Xb(iP/) (J = 0, 1,2), the spin-averaged 



m.e. I(3S\r\lPj) are shown to be small, as it was predicted in a number of studies before. 
However, due to spin-orbit effects the w.f. of the 1 3 Pi state is mixed with the 2P, 3P 



states, for which the m.e. I(3S\r\2P) and I(3S\r\3P) are large and have different signs. 
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Such a mixing is important, although the spin-orbit splittings themselves are not large and 
their typical values are ~ 10—12 MeV. Due to this mixing, a strong cancellation takes 
place in the m.e. I(3S\r\lPi), which gives rise to a suppression of the transition rate for the 
radiative decay T(3S f ) — >■ 7 + %b(l 3 Pi). 

The following partial widths are predicted: Fj(T(3S) — >■ 7 + Xfc(l 3 -Pj)) — 213 eV, 13 eV, 
and 54 eV for J = 2, 1,0, which are in good agreement with the BaBar data, r 2 (exp.) = 
216 ± 25 eV and r (exp.) = 55 ± 10 eV [9]. Also for J = 1 the calculated partial width 
r"i = 13 eV satisfies the upper limit, r\ < 22 eV, obtained in the BaBar experiment. More 
precise measurements of the transition rate for T(3S f ) — > 7 + XbiO-P) could give additional 
restrictions on the spin-orbit effects in radiative decays. 

We predict the following ratio of the partial widths:, r^o = p- = 0.24, which should 
be considered as an important feature of the transition rates where spin-orbit dynamics 
dominates. 
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